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SUMMARY 

A hypothetical turbofan engine simplified simulation with a multivariable 
control and sensor failure detection, Isolation, and accommodation logic 
(HYTESS II) Is presented. The digital program, written In FORTRAN, Is self- 
contained, efficient, realistic, and easily used. Simulated engine dynamics 
were developed from linearized operating point models. However, essential non- 
linear effects are retained. The simulation Is representative of a hypotheti- 
cal, low bypass ratio turbofan engine with an advanced control and failure 
detection logic. Included Is a description of the engine dynamics, the control 
algorithm, and the sensor failure detection logic. Details of the simulation 
Including block diagrams, variable descriptions, common block definitions, sub- 
routine descriptions, and Input requirements are given. Example simulation 
results are also presented. 


INTRODUCTION 

This report Is a user's manual for the hypothetical turbofan engine sim- 
plified simulation with multivariable control and sensor failure detection 
logic (HYTESS II). This simulation builds upon the technology reported In 
reference 1. Essentially, the simulation developed In reference 1 (HYTESS) has 
been modified to Incorporate a control law and sensor failure logic. Addition- 
ally, some Improvements have been Incorporated. In particular some routines 
have been eliminated to simplify the program flow. Also, some Information 
transfer between routines Is now accomplished explicitly using subroutine 
arguments rather than Implicitly using large common blocks. This makes follow- 
ing and understanding program flow much easier. Finally, In HYTESS scaled 
(unitless) variables were used extensively In the simulation along with 
unsealed variables. This Inconsistency often led to confusion. In HYTESS II 
all variables are now unsealed and therefore consistent from one routine to the 
next. Performance scaling, however, Is still a feature of the program. It now 
takes place, once, during program Initialization. 

This digital simulation exists as FORTRAN source code and was designed for 
use on the NASA Lewis Research Center's IBM 3033 AP computer running under the 
TSS/370 operating system. The program Is self-contained and was developed to 
offer those Interested In engine dynamics and controls research an efficient, 
realistic, and easily used engine simulation. 

Typically turbine engine simulations Incorporate detailed nonlinear 
descriptions of both steady-state and dynamic engine operation throughout the 
engine's flight envelope. These detailed nonlinear simulations are very accu- 
rate and realistic and, when Implemented In a digital computer, require 
relatively large amounts of computer storage and computer processing time. 
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This makes these detailed simulations difficult and costly to use. HYTESS II 
was developed as an alternative. It Is structurally simpler than a full non- 
linear engine simulation and therefore has reduced storage and processing 
requirements. HYTESS II retains the essential nonlinear effects Inherent In 
the engine's operation. This Is accomplished by modeling the engine using a 
linear state space formulation, and Incorporating the nonlinear characteristics 
by representing the matrix elements within the linear state space structure as 
nonlinear functions of various engine variables. The compromise Implied In 
this process Is that, although the fidelity of HYTESS II Is maintained for the 
variables considered, It Is very difficult to Identify Individual component 
behavior as In a detailed simulation. Also HYTESS II Is restricted to opera- 
tion In regions about the normal operating line of the engine. The engine 
characteristics simulated by HYTESS II, although hypothetical, are qualita- 
tively similar to those of realistic advance turbofan engines. Typical appli- 
cations for this simulation would Include open-loop engine dynamics studies, 
closed-loop controls analysis, and sensor failure detection performance 
studies. 

This report begins with a description of the engine control and failure 
logic simulated by HYTESS II. Descriptions of the mathematical model of the 
engine and the simulation are given. Finally some results are given. Flow 
charts and variable definitions are also Included. 


SYMBOL LIST 


Acronyms 

ADIA 

HYTESS 

LQR 

MVC 

PI 

SFS 

WSSR 

Variables 

AJ 

ALT 

BLC 

D 

F 

f 


Description 

advanced detection, Isolation, and accommodation 

hypothetical turbofan engine simplified simulation 

linear quadratic regulator 

multivariable control 

proportional -Integral 

sensor failure simulator 

weighted sum of squared residuals 

Description 

nozzle jet area 

altitude 

compressor bleed flow 
feed forward system matrix, 
state system matrix 

engine nonlinear state function vector 
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FGV 

FNMX 

FTIT 

Variables 


G 

9 

H 

JT 

K 

N 

N1 

N2 

PLA 

PT2 

PT4 

PT6 

PO 

S 

SMHC 

SMN 

SVA 

TT2 

TT25 

TT4PLO 

TT45 

TT45LO 

t 


gain system matrix 
fan guide vane angle 
thrust 

fan turbine Inlet temperature 

Description 

control system matrix 

engine nonlinear output function vector 

output system matrix 

log likelihood ratio 

Kalman gain matrix 

number of past residual samples 

fan speed 

compressor speed 

power lever angle 

engine face pressure 

burner pressure 

augmentor pressure 

ambient pressure 

switching matrix 

compressor surge margin 

Mach number 

high compressor stator vane angle 

engine face temperature 

fan discharge temperature 

burner exit slow response temperature 

Inter-turbine temperature 

fan turbine Inlet slow response temperature 

time 
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u 

v 

w 

WFMB 

Variables 

X' 

Y 

6 

c 

0 

X 

$ 

Subscripts 

b 

H 

I 

1 

m 

RP 

S 

SS 

TR 

0 


engine control vector 
value 

residual weighting factor 

main burner fuel flow 

Description 

engine state vector 

engine output vector 

corrected face pressure 

residual vector 

corrected face temperature 

detectlon/lsolatlon threshold 

environmental variables 

Description 

base point 

hard failure 

Isolation 

hypothesis number 

measured 

reference point 

soft failure 

steady-state 

trajectory control 

normal mode 


SYSTEM DESCRIPTION 

The system described by this report consists of a hypothetical turbofan 
engine, sensors, actuators, control system, sensor failure simulator (SFS), and 
sensor failure detection. Isolation, and accommodation (DIA) logic. These com- 
ponents are connected as shown In figure 1. Note that sensors for three 
classes of variables are shown: (1) control variables, (2) environmental 

variables, and (3) engine outputs. Note also that the sensor failure logic 
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applies only to the engine output variable sensors. A description of each of 
these components and Its respective model follows. 


ENGINE DESCRIPTION 

The engine simulated by HYTESS II Is representative of current high tech- 
nology engines and Is shown schematically In figure 2. It Is a low bypass 
ratio, twin-spool, axial-flow turbofan engine, consisting of the following 
components: 

(1) Low-speed fan driven by a turbine (spool 1) 

(2) High-speed compressor driven by a turbine (spool 2) 

(3) Main burner 

(4) Annular fan duct that surrounds the basic gas generator and 
discharges air Into the mixed flow augmentor 

(5) Variable area nozzle 

Variable Inlet guide vanes are used ahead of the fan to Improve Inlet distor- 
tion tolerance and fan efficiency. Variable stators In the high compressor 
Improve starting and high Mach number characteristics. Airflow bleed Is 
extracted at the compressor exit to Improve starting. The exhaust nozzle vari- 
able geometry enables all three nozzle performance parameters (nozzle area, 
expansion ratio, and boattall drag) to be simultaneously near optimum through- 
out the operating range. A list of engine Inputs and outputs Is given In the 
next section. 


Engine Model 

A detailed nonlinear engine model can be written In vector differential 
equation form 


» 

X = f(X,U,S) 

y = g(x,u,$) 


0 ) 


where X Is a state vector, U Is the vector of controls, Y Is the output 
vector, and $ Is a vector of environmental conditions. Detailed nonlinear 
engine relations are represented by the functions f(-) and g(*)- At a 
base point, that Is a steady-state point on the operating line, 


f < vw = 0 


Y b ■ g( VVV 


( 2 ) 


In HYTESS the state space description of the model of equations (1) and (2) Is 
Implemented as 
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( 3 ) 


X s$ = X b - F~ 1 G( U - U b ) 

X = F(X - X $s ) 

Y = Y b + H(X - X b ) f D(U - U b ) 

The subscript b Is used to denote base points. The subscript ss Is used 
to denote the steady-state value of X for a given U. The matrices F, 

F _1 G, H, and D are the typical system matrices. The states, Inputs, and 
outputs were chosen to be typical of those variables used In dynamics and con- 
trols analysis In modern turbofan engines and consist of the following 
variables . 

States : 


X] fan speed (Nl), rpm 
X2 compressor speed (N2), rpm 

X3 burner exit slow response temperature (TT4PL0), K 
X4 fan turbine Inlet slow response temperature (TT45L0), K 

Engine Inputs : 

U-] main burner fuel flow (WFMB), kg/sec 

U2 nozzle jet area (AJ), m 2 

U3 fan guide vane position (FGV), deg 

U4 high compressor variable stator vane angle (SVA), deg 

U5 customer compressor bleed flow (BLC), percent 

Engine outputs : 

Y-] fan speed (Nl), rpm 

Y2 compressor speed (N2), rpm 

Y3 burner pressure (PT4), N/m 2 

Y4 augmentor pressure (PT6), N/m 2 

Y5 fan turbine Inlet temperature (FTIT) , K 

Y 6 thrust (FNMX) , N 

Y7 compressor surge margin (SMHC) 

Operating conditions : 

$1 ambient pressure (PO) 

§2 engine face pressure (PT2), N/m 2 

$3 engine face temperature (TT2), K 

$4 fan discharge temperature (TT25) 

The system matrices were determined In the following manner. Linearized system 
matrices at several base points were found from a representative detailed non- 
linear simulation using perturbatlonal techniques. The elements of each of 
these matrices were regressed upon selected engine variables or elementary 
functions of these variables (elements of Y and $) . As a result, non- 
linear polynomial functions were found that fit the change In these matrix 
elements for the full range of engine power through the flight envelope as 
shown In figure 3. An example of some typical regression polynomials for 
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the system matrices Is given In table 1. Rewriting equation (3) with a more 
explicit functional notation yields 

X ss = x b (Y ’ $) " [F- 1 G](Y,*)[U - U b (Y,#)] \ 

X = F(Y,*)[X - X ss ] ( (4) 

Y = Y b (Y,$) + H( Y ,$) [X - X b (Y,*)] ♦ D(Y,*)[U - U b (Y,*)] \ 

Note, that In the list of engine variables, TT25 was defined as an operating 
condition variable. Although strictly an engine output, TT25 ($ 4 ) Is being 
called an environmental or operating condition variable for three reasons. 

First the only place TT25 Is used In the control logic Is to schedule the SVA. 
Thus, TT25 Is used only In a fashion similar to that for the other environmen- 
tal variables. Second, TT25 Is not covered by the sensor failure logic as are 
all of the listed engine output variables. Finally, In this simulation TT25 
Is modeled as 


TT25 = ^ 0.2308 0.82 j TT2 

Thus, TT25 Is linearly related to TT2, an environmental variable. 


Control Description 

The control law used here Is a multivariable proportional-integral (PI) 
control law. The control used Is a modification of an existing control 
(ref. 2) designed for a high performance turbofan engine. This engine has a 
similar Input-output structure to HYTESS II, and thus, the existing control was 
a logical starting point for the HYTESS II control. A block diagram of the 
HYIESS II control Is given In figure 4. Also shown are the subroutines that 
correspond to the various blocks. The reference point schedules establish the 
desired steady-state performance of the engine and essentially are an Implemen- 
tation of equation 2. The reference point schedules are followed by the tran- 
sition control which limits the rate of change of the commanded variables 
during changes In operating point. The variables, Ujr, serve as feed forward 
controls, while Yjr serves as a command trajectory to be followed by the con- 
trol law. The PI control acts upon differences between Yjr and the engine 
output estimates which are obtained from the DIA logic. The PI gains are 
scheduled as a function of $ and the compressor speed component of Yjr. 
Finally, the maximum and minimum values of the control are limited by the 
engine protection logic to ensure safe engine operation. 


Sensor Failure Logic Description 

Also, Incorporated In the simulation Is logic based upon the principle of 
analytical redundancy (ref. 3) to detect. Isolate, and accommodate (DIA) sensor 
failures. The DIA logic Incorporated In this simulation Is a modification of 
a recently developed algorithm (ref. 4) called the Advanced Detection, Isola- 
tion, and Accommodation (AOIA) algorithm. The ADIA as originally developed 
Incorporates advanced filtering and detection methodologies. The ADIA logic 
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was modified for use with HYTESS II by Incorporating the HYTESS II model In the 
ADIA logic. Thus, In the simulation there are two Identical copies of the 
engine dynamics. The first engine model Is combined with a multivariable con- 
trol and DIA logic to form the simulation. The DIA logic portion of the sim- 
ulation Itself contains the second copy of the engine model which Is used to 
estimate engine outputs. 

The DIA logic consists of four elements: (1) hard failure detection and 

Isolation logic, (2) soft failure detection logic, (3) soft failure Isolation 
logic, and (4) an accommodation filter. The algorithm detects two classes of 
sensor failures, hard and soft. Hard failures are out-of-range or large bias 
errors that occur Instantaneously In the sensed values. Soft failures are 
small bias errors or drift errors that accumulate relatively slowly with time. 
The algorithm Inputs are the measured engine Inputs, U m (t), and the measured 
engine outputs, Y m (t). The algorithm outputs are optimal estimates, Y(t), 
of the engine outputs, Y(t). 

The algorithm has two modes of operation, normal and failure. During 
normal mode operation, l.e., when no sensor failure Is present, the normal mode 
accommodation filter uses all the measured Information to determine Y(t). In 
failure mode operation, one of the five sensors has failed. Simultaneous 
multiple sensor failures are rare events and are not considered. A threefold 
process takes place once the failure has occurred. First the failure Is 
detected. Once a failure Is known to have occurred, the specific faulty sensor 
must be Isolated. Finally, when Isolation has occurred, the failure Is accom- 
modated by recojjf Igurlng the normal mode accommodation filter which generates 
the estimates, Y(t). This threefold procedure takes place for both hard and 
soft failures. 

The normal mode accommodation fljter logic, shown In figure 5, generates 
the estimates of the engine outputs, Y(t). In the Kalman filter equations, the 
matrices F, G, H, and D are typical state space system matrices where X(t) 

Is the 4 by 1 vector of estimates of the engine's state variables and e(t) Is 
the 5 by 1 vector of residuals. The matrix K Is the Kalman gain matrix, and 
S Is a switching matrix. The diagonal elements of S, s^, are either 1 or 
0. All the system matrices as well as the Kalman gain matrix are scheduled as 
a function of operating point to model variations In engine dynamics. 

The hard failure detection and Isolation logic, shown In figure 6, per- 
forms a straightforward threshold check on each sensor residual, c\. Thresh- 
old values are determined from sensor and process noise values as well as sen- 
sor range considerations. If a residual value Is greater than the threshold, 

X^l , hard failure detection and Isolation follow Immediately . 

If a hard failure has not occurred, then a soft failure check Is per- 
formed. The soft failure detection logic, shown In figure 7, first calculates 
an average weighted sum of squared residuals (WSSR). A soft failure Is 
detected when the weighted sum Is greater than a prespecified soft failure 
detection threshold, \$. The number, N, of past residuals summed to obtain 
the average, the weighting factor, W, and the detection threshold, \$, are 
design parameters that are chosen to provide an acceptable tradeoff between 
false alarms and missed detections. 

Once a soft failure has been detected, the soft failure Isolation logic, 
shown In figure 8, Is used to Isolate the failed sensor. Six different 
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Isolation filters generate six different sets of residuals, one for each pos- 
sible failed sensor, and one based on no failed sensors. A log likelihood 
ratio, Is generated for each set of residuals. A test Is then performed 
which determines the most probable set of residuals by finding the maximum 
jtf. When this maximum Is above an Isolation threshold, Xj, the 
faulty sensor Is Isolated. 


Once a hard or soft failure Is detected and Isolated, the accommodation 
filter Is reconfigured by an appropriate change of Its Kalman gain matrix, K, 
to remove the failed sensor from consideration. This Is shown conceptually as 
a modification of the switching matrix S. For example If a failure In sensor 
3 has been Isolated then $33= 0 and s^ = 1, 1 = 1,2, 4, 5. That Is 


S 



0 

1 

0 

0 

0 


0 

0 

0 

0 

0 


0 0 
0 0 
0 0 
1 0 
0 1 


For a soft failure of senior 1, the^accommodatlon filter Is also reinitialized 
to the current value of Y^(t) and X^(t) from the appropriate Isolation 
filter. Reinitialization Is necessary since a significant amount of time may 
have elapsed between failure and Isolation. 


The output of the algorithm, Y, Is directed to the PI control law. Thus, 
although the filter reconfigures due to sensor failures no control mode recon- 
figuration Is required. All DIA logic Is contained within the subroutine FDIA. 
Filter dynamics are calculated by the subroutine FILTER. 


Program Description 

HYTESS II contains all the subroutines necessary to execute the program. 
There are no system library routines required. The program Is written In 
FORTRAN 66, and It Is anticipated that few modifications, If any, will be 
required to execute the program on any system of adequate size that supports 
FORTRAN. Additionally, the program Is entirely compatible with the ANSI 77 
FORTRAN standard except for the use of the NAMELIST feature In the following 
subroutines 6CNTL, INCNTL, MODEL, READIN, SCURVE, and TRCNTL. The NAMELIST 
statements In all the routines, except READIN, could be trivially removed with- 
out changing program operation. Changes to eliminate the use of the NAMELIST 
feature In READIN would be straightforward but substantial. 

The program Itself consists of a main program and 46 subroutines and 5 
block data routines. (See table II for program hierarchy.) There are a total 
of eight levels with a maximum of seven levels of nested subroutines. For 
example at level III the subroutine INLET (called by STDST8) calls four sub- 
routines: ALTABL, PRCMB, HFTA, and TFHA. Subsequently, PRCMB, HFTA, and TFHA 

all call PVAL. Subroutines are listed In the order of their first occurrence 
in the calling program. Several of the subroutines may be called more than 
once by the calling routine. No attempt has been made to show multiple calls 
In table II. In table III a description of the purpose of each subroutine and 
block data Is given. 
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Basic program flow Is shown In figure 9. Essentially, the program Ini- 
tializes, reads Input data, calculates a steady-state point, and If required, 
calculates a user-specified transient. The program accepts as steady-state or 
transient Input either of the following two sets of Inputs: 

Input set 1 (IS1) 

Altitude (Alt) 

Mach number (MN) 

Power level angle (PLA) 

Input set 2 (IS2) 

Altitude (Alt) 

Mach number (MN) 

Fuel flow (WFMB) 

Nozzle area (AJ) 

Fan guide van angle (FGV) 

High compressor stator vane angles (SVA) 

Bleed flow ( BLC) 

Input set 1 Is used to simulate engine response to pilot (PLA) requests. 

Engine response to PLA Is the typical mode of operation for the simulation. 
Input set 2 overrides the engine Inputs calculated by the control logic and 
replaces them with user supplied values. This mode of operation might be use- 
ful for a study of engine dynamics, for example. Also Included In the Input 
Information are program control parameters which control steady-state and 
transient execution of the program. Input requirements are completely speci- 
fied In the next section. This Is followed by descriptions of the steady- 
state and transient options of the main program. 


Input Requirements 

HYTESS II uses the FORTRAN Input mode called NAMELIST to accept values for 
Input parameters. There are six namelists used In subroutine READIN to define 
Input parameters: INPUT, SCALE, MVCIN, INTRAN, PLOT, and FSNS. To help 

Illustrate data entry two examples of program Input which correspond to IS1 and 
IS2 are given In tables IV and V, respectively. Note that although program 
output Is available In both English and metric units, program Input must be 
specified In English units. 


Namelist INPUT 

The namelist INPUT Is used to define the steady-state Input as well as for 
program control. All of the variables used to define a steady-state engine 
condition as well as some option control parameters are contained In namelist 
INPUT. The variable names, their default values, and descriptions are given 
In table VI. For example In table IV the namelist INPUT Is used to (1) Indi- 
cate that a transient Is required (TRAN = 1.0); (2) Indicate that plotting 
variables are to be saved (IPLOT = 1); (3) Indicate that sensor failures are 
to be detected (SFAIL = 1.0); (4) define the Initial conditions for IS1(SALT = 
0.0, SMACH = 0.0, SPLA = 20.0); and (5) indicate that the program output Is to 
be In metric units. 
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Namelist SCALE 


The namelist SCALE Is used to change the program's scale factors. These 
scale factors can be changed to modify engine simulated performance. The vari- 
able names, their default values, and descriptions are given In table VII. The 
default values are selected to yield a performance that Is typical of modern, 
high performance turbofan engines. 


Namelist MVCIN 

The namelist MVCIN Is used to control the modification of the multivaria- 
ble control. The variable names, their default values, and descriptions are 
given In table VIII. 


Namelist INTRAN 

The namelist INTRAN Is used to define the Input necessary for an engine 
transient. This namelist Is only read If TRAN = 1.0 In the namelist INPUT. 

All of the parameters entered via this namelist are entered In the following 
array format 

ARRAY = t 1 ,v 1 ,t 2 ,v 2 t ,v 

where ARRAY represents the respective variable array and t-j ,v-| Is a tlme/- 
value pair. Up to seven pairs may be entered for each array variable, l.e., 
n < 7. The INTRAN namelist parameters are accepted as either step or ramp 
Inputs. Since the particular array variable Is specified at n discrete time 
points, the Intervals between time points need to be further specified. This 
Is accomplished by defining step and ramp Inputs. For a step Input the array 
value In any Interval, say t^ to t^ + -|, Is equal to the preceding time 
point array value v-j . For a ramp Input, the array value In any Interval 
lies on a straight line defined by the two points t^,v^ and t^ + ^,vi 4 .-|. 
Variable names, default values, and descriptions for the transient Input name- 
list INTRAN are given In table IX. For example In table IV namelist INTRAN Is 
used to define the print Interval, the user specified PLA transient, and the 
sampling rate of the control. In this case the PLA Input Is specified as shown 
In figure 10. 


Namelist PLOT 

The namelist PLOT Is used to specify plotted variables. This namelist 
will only be read If IPLOT = 1 In the namelist INPUT. Namelist PLOT contains 
three array variables, IPVAR1, IVPAR2, and IPVAR3 as defined In table X. These 
array variables are used to specify the variables to be stored for plotting. 

A certain variable Is specified for plotting by Including Its associated 
Integer . value (channel number) as defined In appendix A In the variable IPVAR1, 
In appendix B In IPVAR2 and In appendix C In IPVAR3. For example In table IV, 
IPVAR1 contains the channel numbers 1, 4, 22, and 23. Referencing these num- 
bers with appendix A, It can be seen that TIME(l), PLA(4), engine fan 
speed(22), and engine compressor speed(23) are specified for plotting. From 
table IV, IPVAR2 specifies sensed fan and compressor speed for plotting while 
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IPVAR3 specifies the accommodation filter residuals that correspond to fan and 
compressor speed. 


Namelist FSNS 

The namelist FSNS Is used to define various performance parameters (detec- 
tion thresholds, e.g.) for the sensor failure detection logic as well as spec- 
ifying simulated failures. Sensor failures are simulated In the subroutine 
FSENS as 


YFAILED = $F*YSENSED + BIAS +■ NOISE (5) 

Default values (SF = 1.0, BIAS = 0.0, NOISE = 0.0) are for the no failure 
(normal mode) case. The variable NAMP Is used to specify noise amplitude. The 
variables that begin with SN are used to define the scale factors (SF In 
eq. 5), and the variables that begin with BI define the bias failure values. 

A complete description is given In table XI. In the IS1 example of table IV a 
bias failure Is simulated on the PT6 sensor. The failure Is plotted versus 
time In figure 11 . 


Steady-State Program 

After the main program Is Initialized and the Input conditions are spec- 
ified, a steady-state operating point must be established. This Is accom- 
plished by a call to subroutine ST0ST8. Basic program flow through this 
subroutine Is shown In figure 12. First the subroutine INLET Is used to cal- 
culate engine face conditions. For example given ALT and SMN (see appendix A), 
INLET calculates P0, TO, PT2, V0, and ETARAM. This subroutine Is based upon a 
table of altitude, pressure, and temperature data taken from the 1962 Geometric 
Standard. Next, PLA Is determined, either explicitly or Implicitly from the 
user specified Input. Then nominal engine operation Is calculated by RPSCH for 
this value of PLA. Next, the multivariable control and the sensor DIA logic 
Is calculated for the nominal engine response by MVCNT. Next, engine model 
base points (eq. 2) and system matrices are calculated by EM0DEL. Finally, the 
actual engine model response Is simulated by SIMUL using the control Informa- 
tion obtained from MVCNT and the parameters obtained from EM0DEL. A conver- 
gence test Is applied to fan speed. If the fan speed predicted by SIMUL (model 
response) Is equal to the fan speed (nominal response) used In the Initial call 
to RPSCH (within 0.01 percent) then convergence to a steady-state point Is 
achieved. The steady-state values of the engine states, controls, and outputs 
are used as Initial conditions In any subsequent transient calculation. 


Transient Program 

After a steady-state operating point Is established, the main program 
determines If a transient calculation was requested (TRAN = 1.0). If 
requested, the transient Is simulated by the subroutine TRANS. A block diagram 
of the basic program flow through TRANS Is given In figure 13. Here the pro- 
gram flow Is similar to that of STDST8. First model base points and system 
matrices are computed (EM0DEL). Next, engine dynamics are simulated (SIMUL). 
The data Is printed If the print time specified by the user Is equal to the 
simulated time. Next, a call to NUTIME updates the simulated time and updates 
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all the user requested transient Input information that was entered through the 
namelist INTRAN. If simulated time is greater than the stop time there is a 
final data print before the return to the main program. Next, new engine face 
conditions are calculated (INLET) and the control and sensor DIA logic are sim- 
ulated (MVCNT). The program then returns to the call to EMODEL and continues 
the iteration through the program. Each iteration represents an update of the 
Euler integration scheme (in SIMUL) and the iteration continues until the spec- 
ified final time is reached. 


Program Output 

Figure 14 shows a sample of some of the printout for the test case of 
table IV. The program first prints the number of iterations required to reach 
a steady-state point. This output is controlled by the subroutine STDST8. The 
main program output is generated by calls to the subroutine PRINT. For both 
steady-state and transient data subroutine PRINT prints the variables from the 
common ENGOUT (appendix A), common MVCOUT (appendix B), and common DIAOUT 
(appendix C). The ENGOUT variables are labeled ENGINE RESPONSE VARIABLES' with 
four subheadings. The MVCOUT variables are labeled 'CONTROL RESPONSE VARI- 
ABLES' with sixteen subheadings, and the DIAOUT variables are labeled 'DIA 
RESPONSE VARIABLES' with seventeen subheadings. The subroutine PRINT uses a 
10-column format. Each column corresponds to a time at which routine PRINT was 
called. The program also creates an unformatted binary data file written on 
unit 10. This file contains all the transient data that was specified by the 
user for plotting. Figures 15 to 17 are time plots of the variables specified 
for plotting In the example of table IV. The plot of PLA Is already given In 
figure 9. Routines for plotting transient data are not Included In the HYTESS 
II program. Thus, the user must Interface his own plotting capability to the 
program. 


CONCLUDING REMARKS 

A hypothetical turbofan engine simplified simulation with control and 
sensor failure DIA logic Is presented. The program Is suitable for dynamics 
and control analysis. The engine simulation Is structurally simpler than a 
detailed performance digital simulation. However, It does retain the essential 
nonllnearltles of the engine and accurately simulates qualitative engine opera- 
tion. The engine Is modeled using a state space structure. Elements within 
the state matrices are defined by polynomials whose Independent variables are 
functions of engine environment and engine operation. Storage and execution 
time requirements are significantly less than those for a detailed nonlinear 
simulation and are quite reasonable for typical dynamics and control analysis 
studies. The control Is a multivariable proportional plus Integral controller. 
The control represents advanced engine control technology. Also simulated Is 
sensor failure detection, Isolation, and accommodation logic. This DIA logic 
Is based upon state of the art decision making and control theories. Thus, 
HYIESS II represents a realistic, technically advanced, test bed for a variety 
of research objectives within the dynamics and controls arena. 
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APPENDIX A 


ENGINE SIMULATION VARIABLES IN ENGOUT 


The following Is the list of engine variables used In the simulation. 
These variables are all In the COMMON called ENGOUT and are also printed as a 
portion of the hard copy output of the program. 

Channel Variable Units Description 

number 


1 

T 

sec 

Time 

2 

ALT 

m 

Altitude 

3 

SMN 


Mach number 

4 

PLA 

deg 

Power level angle 

5 

PO 

N/M 2 

Ambient pressure 

6 

TO 

K 

Ambient temperature 

7 

DPO 

N/m 2 

Adder to ambient pressure 

8 

DTO 

K 

Adder to ambient temperature 

9 

PT2 

N/m 2 

Engine face pressure 

10 

TT2 

K 

Engine face temperature 

11 

VO 

m/sec 

Airspeed at the Inlet 

12 

ETARAM 


Ram efficiency 

13 

WFMBH 

kg/sec 

Fuel flow 

14 

AFCD 

m 2 

Nozzle area 

15 

FGVPOS 

deg 

Fan guide vane angle 

16 

SVAPOS 

deg 

High compressor stator angle 

17 

BLC 

Percent 

Bleed flow 

18 

SNFAN 

rpm 

Fan physical speed, engine state 

19 

SNCOM 

rpm 

Compressor physical speed, engine state 

20 

TT4PL0 

K 

Burner exit slow response temperature, engine state 

21 

TT45PL0 

K 

Fan turbine Inlet slow response temperature, 
engine state 

22 

SNFM 

rpm 

Fan physical speed, engine output 

23 

SNCM 

rpm 

Compressor physical speed, engine output 

24 

PT4 

N/m 2 

Burner pressure, engine output 

25 

PT6 

N/m 2 

Augmentor pressure, engine output 

26 

TT45 

K 

Fan turbine Inlet, engine output 

27 

FNMX 

N 

Thrust, engine output 

28 

SMHC 


Compressor surge margin, engine output 

29 

XTRA1 


Extra dummy variable 

30 

XTRA2 


Extra dummy variable 
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APPENDIX B 


CONTROL SIMULATION VARIABLES IN MVCOUT 


The following Is a list of control variables used In the simulation. 
These variables are all In the common called MVCOUT and are also printed as a 
portion of the hard copy output of the program. 


Channel 

number 

Variable 

Units 

Description 

1 

PLAMV 

deg 

PLA Input 

2 

ALTMV 

m 

Altitude Input 

3 

SMNMV 


Mach number Input 

4 

TT2MV 

K 

Ambient TT2 

5 

PT2MV 

kPa 

Ambient PT2 

6 

PLASN 

deg 

Sensed PLA 

7 

TT2SN 

K 

Sensed TT2 

8 

PT2SEN 

N/m 2 

Sensed PT2 

9 

SMNSEN 


Sensed Mach number 

10 

PMODE 


Extra variable 

11 

SNFSEN 

rpm 

Sensed fan speed 

12 

SNCSEN 

rpm 

Sensed compressor speed 

13 

TT25SN 

K 

Sensed TT25 

14 

FTITSN 

K 

Sensed FTIT 

15 

PT4SEN 

N/m 2 

Sensed PT4 

16 

PT6MSN 

N/m 2 

Sensed PT6 

17 

RINP1 


Extra variable 

18 

RINP2 


Extra variable 

19 

RINP3 


Extra variable 

20 

RINP4 


Extra variable 

21 

WFMBFB 

kq/sec 

Sensed fuel flow 

22 

AJFB 

m 2 

Sensed nozzle area 

23 

SVAVFB 

deg 

Sensed high compressor stator angle 

24 

FGVVFB 

deg 

Sensed fan guide vane angle 

25 

BLCFB 

percent 

Sensed bleed flow 

26 

PLAEST 

deg 

Estimated PLA 

27 

TT2EST 

K 

Estimated TT2 

28 

PT2EST 

N/m 2 

Estimated PT2 

29 

SMNEST 


Estimated Mach number 

30 

PMDEST 


Extra variable 

31 

SNFEST 

rpm 

Estimated fan speed 

32 

SNCEST 

rpm 

Estimated compressor speed 

33 

T25EST 

K 

Estimated TT25 

34 

FTIEST 

K 

Estimated FTIT 

35 

PT4EST 

N/m 2 

Estimated PT4 

36 

PT6EST 

N/m 2 

Estimated PT6 

37 

REST1 


Extra variables 

38 

REST2 


Extra variables 

39 

REST3 


Extra variables 

40 

REST4 


Extra variables 

41 

WFMBES 

kg/sec 

Estimated fuel flow 

42 

AJEST 

m 2 

Estimated nozzle area 

43 

SVAVES 

deg 

Estimated stator vane angle 
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Channel 

Variable 

Units 

Description 

number 





44 

FGVVES 

deg 

Estimated fan guide vane angle 

45 

BLCEST 

percent 

Estimated bleed 


46 

PLASS 

deg 

Steady state schedule value 

of PLA 

47 

TT2SS 

K 

Steady state schedule value 

of TT2 

48 

PT2SS 

N/m 2 

Steady state schedule value 

of PT2 

49 

SMNSS 


Steady state schedule value 

of SMN 

50 

PMDSS 


Extra variable 


51 

SNFSCH 

rpm 

Steady state schedule value 

of fan speed 

52 

SNCSCH 

rpm 

Steady state schedule value 

of compressor speed 

53 

T25SCH 

K 

Steady state schedule value 

of TT25 

54 

FTISCH 

K 

Steady state schedule value 

of FTIT 

55 

PT4SCH 

N/m 2 

Steady state schedule value 

of PT4 

56 

PT6MSH 

N/m 2 

Steady state schedule value 

of PT6 

57 

TT4P 

K 

Steady state schedule value 

of TT4PL0 

58 

TT45P 

K 

Steady state schedule value 

of TT45PL0 

59 

FNMXSH 

N 

Steady state schedule value 

Of THRUST 

60 

SMHCSH 


Steady state schedule value 

of compressor surge 




margin 


61 

WFMBSH 

kg/sec 

Steady state schedule value 

of fuel flow 

62 

AJSCH 

m 2 

Steady state schedule value 

of nozzle area 

63 

FGVVSH 

deg 

Steady state schedule value 

of fan guide vane angle 

64 

SVAVSH 

deg 

Steady state schedule value 

of stator vane angle 

65 

BLCSH 

percent 

Steady state schedule value 

of bleed 

66 

PLA1R 

deg 

Transition control value of 

PLA 

67 

TT2TR 

K 

Transition control value of 

TT2 

68 

PT2TR 

N/m 2 

Transition control value of 

PT2 

69 

SMNTR 


Transition control value of 

Mach number 

70 

PMDTR 


Extra variable 


71 

SNFTR 

rpm 

Transition control value of 

fan speed 

72 

SNCTR 

rpm 

Transition control value of 

compressor speed 

73 

T25TR 

K 

Transition control value of 

TT25 

74 

FTITR 

K 

Transition control value of 

FTIT 

75 

PT4TR 

N/m 2 

Transition control value of 

PT4 

76 

PT6M1R 

N/m 2 

Transition control value of 

PT6M 

77 

RTR1 


Extra variable 


78 

RTR2 


Extra variable 


79 

RTR3 


Extra variable 


80 

RTR4 


Extra variable 


81 

WFMBTR 

kg/sec 

Transition control value of 

fuel flow 

82 

AJTR 

m 2 

Transition control value of 

nozzle area 

83 

FGVVTR 

deg 

Transition control value of 

fan guide vane angle 

84 

SVAVTR 

deg 

Transition control value of 

stator vane angle 

85 

BLCTR 

percent 

Transition control value of 

bleed 

86 

DELN1 

rpm 

Regulator error, fan speed 


87 

DELN2 

rpm 

Regulator error, compressor 

speed 

88 

DELPT6 

N/m 2 

Regulator error, PT6 


89 

OELTIT 

K 

Regulator error, FTIT 


90 

DELPT4 

N/m 2 

Regulator error, PT4 


91 

DLQWF 

kg/sec 

Regulator control contribution, fuel flow 

92 

OLQAJ 

m 2 

Regulator control contribution, nozzle area 

93 

OLQSVA 

deg 

Regulator control contribution, stator vane angle 

94 

DLQFGV 

deg 

Regulator control contribution, fan guide vanes 
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Channel Variable Units, 
number 


Description 


95 DLQBLC percent Regulator control contribution, bleed 

96 XINFER rpm Integral error, fan speed 

97 XIEPER — Integral error, EPR 

98 XISVER deg Integral error, stator vane 

99 XIFGER deg Integral error, fan guide vanes 

100 XIBLER percent Integral error, bleed 

101 XIFTER K Integral error, FTIT 

102 XIPHER N/m 2 Integral error, high burner pressure 

103 XIPLER N/m 2 Integral error, low burner pressure 

104 XTRAA1 Extra variable 

105 XTRAA2 Extra variable 

106 DITWF kg/sec Euler difference value of WFMBH 

107 0ITAJ m^ Euler difference value of AJCD 

108 DITSVA deg Euler difference value of SVAPOS 

109 DITFGV deg Euler difference value of FGVP0S 

110 DITBLC percent Euler difference value of BLC 

111 XIWFMB kg/sec Integral control contribution of WFMBH 

112 XIAO m2 Integral control contribution of AJCD 

113 XISVAV deg Integral control contribution of SVAPOS 

114 XIFGVV deg Integral control contribution of FGVPOS 

115 XIBLC percent Integral control contribution of BLC 

116 WFMBCL kg/sec Unlimited control output for WFMBH 

117 AJCL m2 Unlimited control output for AJCD 

118 SVAVCL deg Unlimited control output for SVAPOS 

119 FGVVCL deg Unlimited control output for FGVPOS 

120 BLCCL percent Unlimited control output for BLC 

121 WFC0M kg/sec Commanded control output for WFMBH 

122 AFC0M m2 Commanded control output for AJCD 

123 SVAVCM deg Commanded control output for SVAPOS 

124 FGVVCM deg Commanded control output for FGVPOS 

125 BLCCM percent Commanded control output for BLC 

126 WFABSH Unused variable 

127 WFABTR Unused variable 

128 WFABCL Unused variable 

129 WFABCM Unused variable 

130 WFABFB Unused variable 

131 WFABES Unused variable 

132 XSEGLT Unused variable 

133 XL0D Unused variable 

134 DLQWAB Unused variable 

135 XIWFAB Unused variable 

136 XIABER Unused variable 

137 DITWAB Unused variable 

138 ABS5 Unused variable 

139 ABS6 Unused variable 

140 ABS7 Unused variable 

141 XWFFLG Fuel flow limit flag (1 = Max limit, -1 = Min limit 

142 XAJFLG AJ limit Flag (1 = Max limit, -1 = min limit) 

143 XRCFLG SVAPOS limit flag (1 = Max limit, -1 = Min limit) 

144 XCVFLG FGVPOS limit flag (1 = Max limit, -1 = Min limit) 

145 XBLFLG BLC limit flag (1 = Max limit, -1 = Min limit) 

146 XABFLG Unused variable 
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* 


Channel Variable Units Description 

number 

147 XNFFLG Integral control (=1, fan speed control) 

148 XFIFLG Integral control flag (=1, FTIT control) 

149 XPHFLG Integral control flag (=1, Max PT4 control) 

150 XPLFLG Integral control flag (=1, Min PT4 control) 

151 XMTRAN Transient Indication flag 

152 SPR2 Unused variable 

153 SPR3 Unused variable 

154 SPR4 Unused variable 

155 SPR5 Unused variable 

156 SPR6 Unused variable 

157 SPR7 Unused variable 

158 SPR8 Unused variable 

159 SPR9 Unused variable 

160 SPR10 Unused variable 

161 SPR11 Unused variable 

162 SPR12 Unused variable 

163 SPR13 Unused variable 

164 SPR14 Unused variable 

165 SPR15 Unused variable 

166 WASH kg/sec Maximum airflow 

167 WFMIN kg/sec Minimum fuel flow 

168 WFMAX kg/sec Maximum fuel flow 

169 WABMX Unused variable 

170 WABMN Unused variable 

171 FTITMX K Maximum FTIT 

172 PBMIN N/m 2 Minimum PT4 

173 PBMAX N/m 2 Maximum PT4 
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APPENDIX C 


DIA SIMULATION VARIABLES IN DIAOUT 


The following Is a list of DIA variables used In the simulation. These 
variables are all In the common called DIAOUT and are also printed as a 
portion of the hard copy output of the program. 

Channel Variable Description 

number 

1 to 11 SF Scale factor failure values for the eleven sensed variables 

12 to 22 B1 Bias failure values for the eleven sensed variables 

23 to 33 V Noise failure values for the eleven sensed variables 

34 to 38 Z The vector of measurements used as the accommodation 

filter Input 

39 to 43 ZUL The vector of unlagged (no sensor dynamics) estimates of Z 

44 to 48 ZH The vector of estimates of Z (with sensor dynamics) 

49 to 53 GAMDHO The vector of accommodation filter residuals 

54 DTHRSH Unused 

55 DFLAGH Hard detection flag 

56 KI Normalized and squared sum of residuals 

57 ETA Weighted sum (over time) of squared residuals 

58 DTHRSS Unused 

59 DFLAGS Soft detection flag 

60 to 64 ZIO Output vector of hypothesis 0 

65 to 69 GAMIHO Residual vector of hypothesis 0 

70 to 74 ZI1 Output vector of hypothesis 1 

75 to 79 GAMIH1 Residual vector of hypothesis 1 

80 to 84 ZI2 Output vector of hypothesis 2 

85 to 89 GAMIH2 Residual vector of hypothesis 2 

90 to 94 ZI3 Output vector of hypothesis 3 

95 to 99 GAMIH3 Residual vector of hypothesis 3 

100 to 104 ZI4 Output vector of hypothesis 4 

105 to 109 GAMIH4 Residual vector of hypothesis 4 

110 to 114 ZI5 Output vector of hypothesis 5 

115 to 119 GAMIH5 Residual vector of hypothesis 5 

120 HI00 Likelihood value for hypothesis 0 

121 to 125 HI Likelihood value for hypothesis 1 to 5 

126 to 130 H0MHI Log likelihood ratios 

131 ITHRSH Unused 

132 to 136 ICHAN Isolated channel flags 

137 IS0LT Isolation flag 

138 ACCOM Accommodation flag 

139 to 145 ZACC0M Accommodation filter outputs 

146 to 200 D0UT Extra variables 
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TABLE I. - TYPICAL REGRESSION POLYNOMIALS FOR THE ELEMENTS OF 

THE SYSTEM MATRICES 


System 

matrix 

Element 

Polynomial 3 

F 

(1,1) 

0.0968 0.0019 PT6 2 „ 

- * - - 2 ' 463 

F _1 G 

(4,1) 

0.000933 PT6 0.97xl0~ 9 N1 PT6 2 n n „ rn . 

c — • «"■ — U.UJbUb 

* V? i 2 

H 

(5,3) 

0.03110 * 0.5486x1 0 -4 - 0.3612 

D 

(5,1) 

-0.03540 + 31 p j 256 + 0.04914 


a S = Pl/14.696; e = Tl/518.67 
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TABLE II. - HYTESS II ROUTING HIERARCHY 


Level 


I 

II 

III 

IV 

V 

VI 

VII 

VIII 

MAIN 

SETUP 

MSET 

READIN 

SCURVE 

UNBAR 







GCNTL 

UNBAR 







TRCNTL 

MODEL 

UNBAR 

UNBAR 






INCNTL 

UNBAR 






ST0ST8 

INLET 

ALTABL 

PRCMB 

PVAL 







HFTA 

PVAL 







TFHA 

PVAL 






SNFMAP 

UNBAR 







SN1SCH 

SCURVE 

UNBAR 






N2TABL 

SCURVE 

UNBAR 






RPSCH 

N2TABL 

SCURVE 

UNBAR 






SCURVE 

SPRINT 

FGVCAL 

UNBAR 






MVCNT 

SENSOR 

XLAGX 







FSENS 

GRAND 

RANDU 






EMODEL 

SNFMAP 

UNBAR 







RPSCH 

N2TABL 

SCURVE 

UNBAR 






SCURVE 

SPRINT 

FGVCAL 

UNBAR 





FDIA 

FILTER 

SUB 

MUL 

ADD 

DMINV 

SCA 

XLAGX 







NORMAL 

USSR 

XLAGX 

MUL 






RPSCH 

N2TABL 

SCURVE 

UNBAR 






SCURVE 

SPRINT 

FGVCAL 

UNBAR 






TRCNTL 

MODEL 

UNBAR 

UNBAR 






GCNTL 

UNBAR 







LQR 

INCNTL 

UNBAR 







EPRTKT 

ACTUAT 

SCURVE 

UNBAR 





EMOOEL 

SNFMAP 

UNBAR 







RPSCH 

N2TABL 

SCURVE 

UNBAR 






SCURVE 

UNBAR 

















TABLE II. - Concluded. 


Level 


III 

IV 

V 

VI 

VII 

VIII 

SIMUL 

SUB 

MUL 

SCA 

ADD 

SPRINT 

FGVCAL 





PRINT 




NUTIME 




INLET 

ALTABL 




PRCMB 

PVAL 



HFTA 

PVAL 



TFHA 

PVAL 


MVCNT 

SENSOR 

XLAGX 



FSENS 

GRAND 

RANDU 


EMODEL 

SNFMAP 

UNBAR 



RPSCH 

N2TABL 




SCURVE 




SPRINT 




FGVCAL 


FOIA 

FILTER 

SUB 




MUL 




ADD 




DMINV 




SCA 




XLAGX 



NORMAL 

MUL 



WSSR 




XLAGX 



RPSCH 

N2TABL 

SCURVE 



SCURVE 

UNBAR 



SPRINT 




FGVCAL 



TRCNTL 

MODEL 

UNBAR 


GCNTL 

UNBAR 



LQR 




INCNTL 

UNBAR 



EPRTKT 

SCURVE 

UNBAR 


ACTUAT 



EHODEL 

SNFMAP 

UNBAR 



RPSCH 

N2TABL 

SCURVE 



SCURVE 

UNBAR 



SPRINT 




FGVCAL 


SIMUL 

SUB 




MUL 




SCA 




ADD 




SCURVE 

UNBAR 


UNBAR 


UNBAR 


UNBAR 













TABLE III. - HYTESS II ROUTINE DESCRIPTIONS 

ACTCRV - Block Data, contains data for stator vane angle schedule 
ACTUAT - Simulates actuators 
ADD - General matrix addition 

ALTABL - Calculates temperature and pressure corrections at various altitudes 
DIAPNT - Block Data, contains names of DIA variables to be printed 
DMINV - General matrix Inversion 

EHODEL - Calculates engine model matrices and basepolnts and Kalman filter 
gain matrix 

ENGPRN - Block Data, contains names of engine variables to be printed 

EPRTKT - MVC engine protection logic 

FDIA - DIA analytical redundancy algorithm 

FGVCAL - Calculates fan guide vane angle scheduled value 

FILTER - DIA Kalman filter calculation 

FSENS - Simulates sensor failures 

GCNTL - Calculates proportional and Integral gains 

GRAND - Converts uniformly distributed random numbers to Gaussian distributed 
random numbers 

HFTA - Calculates enthalpy as a function of temperature 
INCNTL - Integral control logic 

INLET - Solves for engine Inlet conditions from ambient conditions 
LQR - Linear Quadratic Regulator (proportional) control logic. 

MAIN - Main program for HYTESS II 

METRIC - Converts from English to metric units or vice versa 
MODEL - Used by the transition control to generate trajectory value of states 
and outputs 

MSET - Initializes program control parameters 

MUL - General matrix multiply 

MVCGAN - Block Data, contains data for PI gain matrices 

MVCNT - Control system structure Including sensors, DIA logic, multivariable 
PI logic, and actuators 

MVCPRN - Block Data, contains names of control variables to be printed 
NORMAL - Computes normalized dot product of a vector 
NUTIME - Finds print Interval and values of ramp functions 
N2TABL - Calculates N2 as a function of PLA and TT2, or PLA as a function of 
N2 and TT2 

PRCMB - Calculates specific heat as a function of temperature 

PRINT - Formats and prints output 

PVAL - Evaluates a given polynomial 

RANDU - Generates uniformly distributed random numbers 

READIN - Reads Input data 

RPSCH - Calculates steady-state operating point values for state, output, 
and control variables 

SCA - Multiply elements of a matrix by a constant 

SCURVE - Contains data (curves) for steady-state operating points 

SENSOR - Simulates sensors 

SETUP - Initialization routine 

SIMUL - Updates engine state space model equations using Euler Integration 
SNFMAP - Finds PLA as a function of N1 , TT2, and SMN 
SN1SCH - Finds N1 as a function of TT2, SMN, and PLA 
SPRINT - Table lookup routine 



TABLE III. - Concluded. 


ST0ST8 

SUB 

TFHA 

TRANS 

TRCNTL 

UNBAR 

WSSR 

XLAGX 


Controls execution during steady-state convergence 
General matrix subtraction 

Calculates temperature as a function of enthalpy 
Controls execution during a transient. 

Transition control logic 
Table lookup routine 

DIA weighted sum of squared residuals calculation 
First order lag 


TABLE IV. - IS1 EXAMPLE INPUT TEST CASE 


&INPUT 

TRAN=1 .0, IPL0T= 1,SFAIL=1.0, 

SALT =0.0, SMACH=0.0,SPLA=20.0 
IMTRC= 1, 

SEND 

&INTRAN 

PNTBLK = 0.0, 0.1, 10.0, .1, 

PLABLK = 0.0,20.0, .5, 20.0, 1.0,83.0, 

SMPBLK = 0.0,0.02, 0.04,0.02, 10.0, .02, 

SEND 
&PL0T 

IPVAR1 = 1, 4,22,23 

IPVAR2 = 11,12 

IPVAR3 = 49,50 

&END 
&FSNS 

Bill BK = 0.0, 0.0, 2. 5, 0.0, 7.502,250., 

SEND 


10.0,83.0, 


TABLE V. - IS2 EXAMPLE INPUT TEST CASE 


&INPUT 

SMACH=0 . ,SALT=0. ,SPLA=52 . ,TRAN=1 . , 

SFGVV=-25. ,SSVAV=6.0, 

SEND 

&INTRAN 

PNTBLK=0. , .1 ,10. , .1 , 

FGVBLK=0. 0,-25. 0,0.1 ,-25.0, .2,-22.5,10. ,-22.5, 
SVABLK=0. 0,6. 0,5. 0,6. 0,5.1 ,5.4,10. ,5.4, 

&END 






TABLE VI. - DESCRIPTION OF STEADY-STATE NAMELIST INPUT 


Variable 

name 

Default 

Description 

TRAN 


If TRAN=1 then transient run desired 

I PLOT 


If I PLOT =1 then plotting desired 

ISCL 


If ISCL=1 then enter new engine scale factors 

SFAIL 


If SFAI L=1 then Include sensor DIA logic 

MFLAG 


If MFLAG=1 then enter new control logic parameters 

SALT 

-999 

Altitude, m 

SMACH 


Mach number 

SPLA 


Power level angle, deg 

STAM 


Ambient temperature, K 

SPAM 


Ambient pressure, N/m 2 

SDTAM 


Adder to ambient temperature, K 

SDPAM 


Adder to ambient pressure, N/m 2 

SPT2 


Engine face total pressure, N/m 2 

STT2 


Engine face total temperature, K 

SWF 


Fuel flow, kg/sec 

SAJ 


Nozzle jet area, m 2 

SFGVV 


Fan guide vane angle, deg 

SSVAV 


Compressor stator vane angle, deg 

SBLC 

T 

Bleed flow, percent 

IMTRC 

0 

If IMTRC=1 then convert output from English to 
metric units 


TABLE VII. - DESCRIPTION OF NAMELIST SCALE 


Variable 

name 

Default value 

Variable description 

1 

10 000 rpm 

Scale factor for engine state l(SNFAN) 

m r 4sQlwm 

15 000 rpm 

Scale factor for engine state 2(SNC0M) 


1 600 °F 

Scale factor for engine state 3(TT4PLQ) 

SC L( 4 ) 

1 600 °F 

Scale factor for engine state 4(TT45L0) 

SCL( 5) 

10 000 rpm 

Scale factor for engine output l(SNFN) 

SCL( 6) 

1 5 000 rpm 

Scale factor for engine output 2(SNCM) 

SCL(7) 

550 psl 

Scale factor for engine output 3(PT4) 

SCL( 8) 

130 psl 

Scale factor for engine output 4(PT6) 

SCL(9) 

1 600 °F 

Scale factor for engine output 5(TT45) 

SCL(IO) 

25 000 lb 

Scale factor for engine output 6(FNMX) 













TABLE VIII. - DESCRIPTION OF NAMELIST MVCIN 


Variable 

name 

Default 

Description 

. 

XMVC 

Defined In 

Equlvalenced to the 81 variables In the 


routine MSET 

common MVCCON 

SSCH 

0.0 

If SSCH^O.O then read new steady-state 
schedule data 

GAININ 

0.0 

If GAINING. 0 then read new PI control gains 

TCNIN 

0.0 

If TCNIN^O.O then read new transition 
control rate data 

XINTIN 

0.0 

If XINTIN^O . 0 then read Integral control 
dead zone data 


TABLE IX. - DESCRIPTION OF TRANSIENT NAMELIST INTRAN 


Variable 
array name 

Default 

Description 

Step Input 

PNTBLK 

SMPBLK 

-999 

-999 

Print Interval 

Control sampling Interval 

Ramp Inputs 

PLABLK 

ALTBLK 

XMNBLK 

WFBLK 

AJBLK 

FGVBLK 

SVABLK 

BLCBLK 

1 

PLA ramp Input 

Altitude ramp Input 

Mach number ramp Input 

Fuel flow ramp Input 

Nozzle area ramp Input 

Fan guide vane ramp Input 

Compressor stator vane angle ramp Input 

Bleed flow ramp Input 


TABLE X. - DESCRIPTION OF NAMELIST PLOT 


Variable 

name 

Default 

Variable description 

IPVAR1 

-1 

Integer array that defines those variables In the common 
ENGOUT that are to be saved for plotting. See appendix A 

IPVAR2 

-1 

Integer array that defines those variables In the common 
MVCOUT that are to be saved for plotting. See appendix B. 

1PVAR3 

-1 

Integer array that defines those variables In the common 
DIAOUT that are to be saved for plotting. See appendix C. 
























TABLE XI. - DESCRIPTION OF NAMELIST FSNS 


Variable Default Variable description 

name 

GMEANO(l) Noise mean for DIA variable 1 

(2) Noise mean for DIA variable 2 

(3) Noise mean for DIA variable 3 

(4) Noise mean for DIA variable 4 

(5) ▼ Noise mean for DIA variable 5 

GSD0(1) 350 Noise standard deviation for DIA variable 1 

(2) 400 Noise standard deviation for DIA variable 2 

(3) 30 Noise standard deviation for DIA variable 3 

(4) 5 Noise standard deviation for DIA variable 4 

(5) 250 Noise standard deviation for DIA variable 5 

NHIST 3 The number of past values used In the calculation of 

the WSSR statistic 

THRH 0.2 Heal threshold factor 

THRDI 2.0 Hard failure detectlon/lsolatlon threshold factor 

THRD 1.6 Soft failure detection threshold 

THRI 3.0 Soft failure Isolation threshold 

TAU .1 Time constant used for smoothing hypothesis log 

likelihoods 

FSPRT 0 If FSPRT=0.0 then write the failure detection data 

and failure scenario data 

NAMP(l) 0 Noise amplification factor for PLASN 

(2) Noise amplification factor for TT2SN 

(3) Noise amplification factor for PT2SN 

(4) Noise amplification factor for SMNSN 

(5) Noise amplification factor for EXTRA 

(6) Noise amplification factor for SNFSEN 

(7) Noise amplification factor for SNCSEN 

(8) Noise amplification factor for TT25SN 

(9) Noise amplification factor for FTITSN 

(10) Noise amplification factor for PT4SEN 

(11) Noise amplification factor for PT6MSN 

SN1BLK -999 Data for simulated step change In PLASN scale factor 

SN2BLK Data for simulated step change In TT2SN scale factor 

SN3BLK Data for simulated step change In PT2SN scale factor 

SN4BLK Data for simulated step change In SNMSN scale factor 

SN5BLK Data for simulated step change In EXTRA scale factor 

SN6BLK Data for simulated step change In SNFSEN scale factor 

SN7BLK Data for simulated step change In SNCSEN scale factor 

SN8BLK Data for simulated step change In TT25SN scale factor 

SN9BLK Data for simulated step change In FTITSN scale factor 

SN10BLK Data for simulated step change In PT4SEN scale factor 

SN11BLK * Data for simulated step change In PT6MSN scale factor 

BI1BLK -999 Data for simulated ramp change In PLASN bias factor 

BI2BLK Data for simulated ramp change In TT2SN bias factor 

BI3BLK Data for simulated ramp change In PT2SN bias factor 

BI4BLK Data for simulated ramp change In SMNSN bias factor 

BI5BLK Data for simulated ramp change In EXTRA bias factor 

BI6BLK u Data for simulated ramp change In SNPSEN bias factor 

BI7BLK ’ Data for simulated ramp change In SNCSEN bias factor 







TABLE XI. - Concluded. 


Variable 

name 

Default 

Variable description 

BI8BLK 
BI9BLK 
Bill BLK 
Bill BLK 

-999 

f 

Data for simulated ramp change In TT25SN bias factor 
Data for simulated ramp change In FTITSN bias factor 
Data for simulated ramp change In PT4SEN bias factor 
Data for simulated ramp change In PT6MSN bias factor 








Station 


ALT SMN 



PLA 

Figure 1. - HYTESS II system block diagram. 
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Figure 2. - Schematic representation of a hypothetical turbofan engine. 











Figure 3. - Engine flight envelope with engine operating 
poi nts. 
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Figure 4. - HYTESS II control block diagram. 
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Figure 5. - Sensor DIA normal mode accommodation filter. 
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Figure 6. - DIA hard failure detection Figure 7. - DIA soft failure detection 

and isolation logic. logic. 
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Figure 8. - DIA soft failure isolation logic. 













Figure 9. - Basic program flow through MAIN. 












Figure 11. - PT6 


MAIN 




Figure 12. - Basic program flow through STDST. 






























Figure 13. - Basic program flow through TRANS. 










HYPOTHETICAL TURBOFAN CONTROLLER 


AAA ENGINE RESPONSE VARIABLES aaa 


1 TIME 

0.00000 

0.10000E 

00 0.20000 

0.30000 0.60000 0.50000 

***AMBIENT CONDITIONS* 

0.60000 

0.70000 

0.80000 

0.90000 

2 ALT 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

3 SMN 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

6 PL A 

20.000 

20.000 

20.000 

20.000 

20.000 

20.000 

32.600 

65.200 

57.800 

70.600 

5 PO 

16.696 

16.696 

16.696 

16.696 

16.696 

16.696 

16.696 

16.696 

16.696 

16.696 

6 TO 

518.69 

518.69 

518.69 

518.69 

518.69 

518.69 

518.69 

518.69 

518.69 

518.69 

7 DPO 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

8 OTO 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

9 PT2 

16.696 

16.696 

16.696 

16.696 

16.696 

16.696 

16.696 

16.696 

16.696 

16.696 

10 TT2 

518.67 

518.67 

518.67 

518.67 

518.67 

518.67 

518.67 

518.67 

518.67 

518.67 

11 VO 

0.00000 

0.00000 

0.00000 

0.00000 

0.00 000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

12 ETARAM 

1.0000 

1.0000 

1.0000 

1.0000 

***ENGINE 

1.0000 

INPUTS*** 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

13 UFMBH 

1133.1 

1133.1 

1133.3 

1133.6 

1133.8 

1136.0 

1227.9 

1358.6 

1669.6 

1566.7 

16 A JCD 

3.0000 

2.9500 

2.9500 

2.9500 

2.9500 

2.9500 

3.0062 

3.0086 

3.0058 

2.9995 

15 FGVP05 

-25 .000 

-25.000 

-25.000 

-25.000 

-25.000 

-25.000 

-26.857 

-26.711 

-26.562 

-26.652 

16 SVAPOS 

-39.033 

-39.033 

-39.033 

-39.033 

-39. D29 

-39.018 

-38.821 

-38.608 

-38.281 

-37.768 

17 BIC 

0.00000 

0.00000 

0.00000 

0.00000 

***ENGINE 

0.00000 

STATES*** 

0.00000 

0.00000 

0.00000 

0 .0000 0 

0.00000 

18 SNFAN 

3560.9 

3536.8 

3536.1 

3531.9 

3530.1 

3528.6 

3536.2 

3557.1 

3592.7 

3662.6 

19 SNCOJ1 

8695.6 

8693.8 

8692.5 

8691.6 

8690.6 

8689.6 

8702.2 

8763.5 

8809.1 

8892.1 

20 TT6PIO 

87.559 

87.506 

87.667 

87.639 

87.618 

87.606 

87.813 

89.101 

91.163 

93.731 

21 TT65L0 

61.513 

61.672 

61.650 

61.639 

*x*ENGINE 

61.637 

OUTPUTS*** 

61.639 

61.998 

63.712 

66.290 

69:368 

22 SNFAN 

3560.9 

3536.8 

3536.1 

3531.9 

3530.1 

3528.6 

3536.2 

3557.1 

3592.7 

3662.6 

23 SNCOfl 

8695.6 

8693.8 

8692.5 

8691.6 

8690.6 

8689.6 

8702.2 

8763.5 

8809.1 

8892.1 

26 PT6 

65.698 

65.670 

65.668 

65.629 

65.616 

65.603 

66.663 

68.769 

71.097 

73.666 

25 PT6 

13.516 

13.617 

13.613 

13.610 

13.608 

13.606 

13.520 

13.576 

13.656 

13.758 

26 TT 65 

976.63 

975.06 

975.51 

975.96 

976.29 

976.55 

1015.6 

1065.2 

1100.9 

1125.1 

27 FHMX 

1363.9 

1381.2 

1380.2 

1379.6 

1378.7 

1378.2 

1373.3 

1399.1 

1630.1 

1665.9 

28 SMHC 

0.77921 

0.77891 

0.77861 

0.77836 

0.77809 

0.77785 

0.76993 

0.76686 

0.76730 

0.77696 


Figure 14. - Sample of program printout for test case of table IV. 
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Figure 15. - Fan and compressor speed versus time 
(IPVAR1). 



14 000 



Time, sec 


Figure 16. - Sensed fan and compressor speed versus 
time (IPVAR2I. 
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